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Abstract 

Based on a symmetry analysis of the microscopic Hubbard and t-J 
models, a systematic low-energy effective field theory is constructed for 
hole-doped antiferromagnets on the honeycomb lattice. In the antiferro- 
magnetic phase, doped holes are massive due to the spontaneous break- 
down of the SU (2) s symmetry, just as nucleons in QCD pick up their mass 
from spontaneous chiral symmetry breaking. In the broken phase the ef- 
fective action contains a single-derivative term, similar to the Shraiman- 
Siggia term in the square lattice case. Interestingly, an accidental continu- 
ous spatial rotation symmetry arises at leading order. As an application of 
the effective field theory we consider one-magnon exchange between two 
holes and the formation of two-hole bound states. As an unambiguous 
prediction of the effective theory, the wave function for the ground state 
of two holes bound by magnon exchange exhibits /-wave symmetry. 
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1 Introduction 



The physics of correlated electron systems is strongly influenced by the geometry of 
the underlying crystal lattice. For example, at weak coupling the half-filled Hubbard 
model on the honeycomb lattice is a semi-metal with massless fermion excitations re- 
siding in two Dirac cones. This situation is realized in graphene. At stronger coupling 
the SU(2) S symmetry breaks spontaneously and the system becomes an antiferro- 
magnet, which may be realized in the dehydrated version of Na2CoC>2 x yH 2 0. On 
a square lattice, on the other hand, due to Fermi surface nesting, the system is an 
antiferromagnet even at arbitrarily weak coupling. Upon doping, antiferromagnets on 
both the square and the honeycomb lattice may become high-temperature supercon- 
ductors. Recenty, a spin-liquid phase was identified in numerical simulations between 
the free fermion graphene phase and the strongly correlated antiferromagnetic phase 

Q. 

The low-energy physics of undoped antiferromagnets on a bipartite lattice is de- 
scribed by an 0(3)-symmetric non-linear a-model l2j, whose systematic treatment is 
realized in magnon chiral perturbation theory [3J-I8]. The effective theory for holes 
doped into an antiferromagnet on the square lattice was pioneered by Shraiman and 
Siggia jof. In particular, these authors found an important term in the magnon- hole 
action known as the Shraiman- Siggia term. Based on the microscopic t-J model, in- 



teresting results on magnon-mediated forces between holes were obtained in [10(| and 



spiral phases were studied in |12| . In analogy to baryon chiral perturbation theory 



for QCD [13Hl7l|. a systematic low-energy effective field theory for magnons and holes 



was constructed in [18|, ll9(. This theory has been used in a detailed analysis of two- 



hole states bound by one- magnon exchange 19|,|20| as well as of spiral phases |2l|. The 



systematic effective field theory investigations have also been extended to electron- 
doped antiferromagnets j22[. In that case, no Shraiman-Siggia-type term (with just 
a single spatial derivative) exists. Hence at low energies magnon-electron couplings 
are weaker than magnon-hole couplings. As a consequence, in contrast to hole-doped 
systems, in electron-doped systems there are no spiral phases with a helical structure 
in the staggered magnetization 22 . 



In this paper we construct a systematic low-energy effective field theory for hole- 
doped antiferromagnets on the honeycomb lattice. In the antiferromagnetic phase, 
the SU(2) S spin symmetry is spontaneously broken and the fermions pick up a mass. 
This is analogous to QCD where protons and neutrons pick up their masses due to 
spontaneous chiral symmetry breaking. Our analysis shows that the effective theory 
on the honeycomb lattice contains a term similar to the Shraiman-Siggia term in the 
square lattice case \^ , which supports spiral phases. Remarkably, the leading terms of 
the effective theory have an accidental continuous rotation symmetry which is reduced 
to the discrete 60 degrees rotation symmetry O of the microscopic honeycomb lattice 
only by the higher-order terms. While spiral phases in hole-doped antiferromagnets 



on the honeycomb lattice were explored in [23|, here — as a further application of 
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the effective field theory method — we derive the one-magnon exchange potentials 
between two holes and study the formation of two-hole bound states, which will turn 
out to have /-wave symmetry. 

The rest of the paper is organized as follows. Section [2] contains a symmetry 
analysis of the underlying Hubbard and t-J model. Section [3] is devoted to the low- 
energy effective theory for magnons — in particular, a non-linear realization of the 
SU(2) S spin symmetry is constructed. Based on the microscopic t-J model, in Section 
HI we include the holes in the effective field theory framework. In Section 0, one- 
magnon exchange potentials between two holes are derived and the resulting two-hole 
bound states are investigated in Section El Finally, section [7] contains our conclusions. 



2 Microscopic Theory 

We assume that the Hubbard model and the t- J-model are reliable models to describe 
doped quantum antiferromagnets, and therefore are valid as concrete microscopic 
models for the low-energy effective field theory for magnons and holes. Due to the 
fact that the effective Lagrangian to be constructed must inherit all symmetries of 
the underlying microscopic systems, a careful symmetry analysis of these microscopic 
models is presented in this section. 



2.1 Symmetries of the Honeycomb Lattice 

The honeycomb lattice is not a Bravais lattice — rather, it consists of two triangular 
Bravais sublattices A and B, as depicted in Figure [TJ The primitive vectors that 
generate the triangular sublattices in coordinate space are given by 

a x = V3a f"^>^J > fl 2 = VSa(0,l), (2.1) 

where a is the lattice spacing between two neighboring sites. The two basis vectors 
bi and b 2 that span the reciprocal lattice obey 

dibj = 2n5ij, (2.2) 

and are given by 

The vectors bi and b 2 generate the hexagonal first Brillouin zone of the triangular lat- 
tice. Since the honeycomb lattice consists of two triangular sublattices, its momentum 
space is doubly-covered. 
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Figure 1: Bipartite non-Bravais honeycomb lattice consisting of two triangular Bravais 
sublattices. The translation vectors are a± and a<i- 

The honeycomb lattice exhibits a number of discrete symmetries. Translations by 
the vectors Oj are denoted by -Dj. Counter-clockwise rotations by 60 degrees around 
the center of a hexagon are denoted by O, and reflections at the Xi-axis going through 
the center of the hexagon are denoted by R. Translations by other distance vectors, 
rotations by other angles or around other centers, and reflections with respect to other 
axes can be obtained as combinations of the elementary symmetry operations D ± , D 2 , 
O, and R. 



2.2 Symmetries of the Hubbard Model 

Let c xs denote the operator which creates a fermion with spin s G {t, i} on a lattice site 
x = (xi,X2)- The corresponding annihilation operator is c xs . These fermion operators 
obey the canonical anticommutation relations 

{ C L> c vA = $xyS ss >, {c xs , c ys ,} = {4 S , c\ s ,} = 0. (2.4) 
The second quantized Hubbard Hamiltonian is defined by 
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where (x, y) indicates summation over nearest neighbors, t is the hopping parameter, 
and the parameter U > fixes the strength of the Coulomb repulsion between two 
fermions located on the same lattice site. The parameter // denotes the chemical 
potential. 

The fermion creation and annihilation operators can be used to define the following 
SU(2) S Pauli spinors 

4 =(44), *=(£)■ 
In terms of these operators, the Hubbard model can be reformulated as 

H = ~t J>t c, + c\c x ) + | Cx _ 1)2 _ ^ Y+c\c x - 1). (2.7) 

The parameter /z = //— ^ controls doping where the fermions are counted with respect 
to half-filling. 

Since all terms in the effective Lagrangian must be invariant under all symmetries 
of the Hubbard model, a careful symmetry analysis of Eq. (12. 7p is needed. Let us divide 
the symmetries of the Hubbard model into two categories: Continuous symmetries 
(SU(2) S , U(1)q fermion number and its non-Abelian extension SU(2)q), which are 
internal symmetries of Eq. fl2.7p . and discrete symmetries (Di, O and R), which are 
symmetry transformations of the underlying honeycomb lattice. There is also time 
reversal which is implemented by an anti-unitary operator T. This symmetry will be 
discussed further in the effective field theory framework. 

In order to construct the appropriate unitary transformation representing a global 
SU(2) S spin rotation, we first define the total SU(2) S spin operator by 

£=£& = £ c I§*- ( 2 - 8 ) 

X X 

The spin symmetry is implemented by the unitary operator 

V = exp(iff- S), (2.9) 

which acts on c x as 

c' x = V^c x V = exp(tff-^)c x = gc x , g G SU (2) s . (2.10) 

The total spin is conserved and the Hubbard Hamiltonian is invariant under global 
SU(2) S spin rotations. This symmetry, however, is spontaneously broken: the corre- 
sponding order parameter is the staggered magnetization vector 

m s = £(-ir4 (2.H) 
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which takes a non-zero expectation value in the ground state of the antiferro magnet. 
We define (— l) x = 1 for all x G A and (— l) x = —1 for all x G B, where A and B are 
the two triangular sublattices of the honeycomb lattice. 

The unitary transformation of the U(1)q symmetry involves the charge operator 

q = J2®* = Efe - 1) = 5>V*t + - 1). ( 2 - 12 ) 

XX X 

which counts the fermion number with respect to half-filling. The corresponding 
unitary operator is given by 

W = exp(iuQ), (2.13) 
and the fermion operators transform as 

Q Cx = W^c x W = exp{iu)c x , exp{iu) G U(l) Q . (2.14) 

Charge or fermion number are conserved due to [H, Q] = 0. 

The Hubbard model shows invariance under shifts along the two primitive lattice 
vectors a\ and a 2 . These transformations are generated by the unitary operators Z?j, 
which act on the spinor c x as 

Dl c x = D\c x Di = c x + ai . (2.15) 

By applying Eq. (l2.15j) on the Hubbard Hamiltonian of Eq. (12.71) and redefining the sum 
over lattice sites x, one can see that indeed [H, Dj\ = 0. Since the shift symmetry 
maps A — y A and B — > B, this transformation does not affect the order parameter 
M s - 

A spatial rotation by 60 degrees leaves Eq. fl2.7p invariant. Since spin-orbit cou- 
pling is neglected in the Hubbard model, spin decouples from the spatial motion and 
becomes an internal quantum number. The rotation symmetry is implemented by the 
use of a unitary operator O, which acts on the fermion operators as 

°c x = 0^c x O = c Qx . (2.16) 

Rotation symmetry on the honeycomb lattice is spontaneously broken because O 
exchanges the two sublattices A -h- B and therefore the staggered magnetization M s 
gets flipped. This is, however, just the same as redefining the sign of (— l) x and does 
therefore not change the physics. In the construction of the effective field theory 
for magnons and holes, it will turn out to be useful to also consider the combined 
symmetry O' consisting of a spatial rotation O and a global SU(2) S spin rotation 
g = ia 2 . O' transforms Cx clS 

°'c x = 0'^c x O' = (ia 2 ) °c x = (ta 2 )c 0x . (2.17) 

The specific SU(2) S element g = ia 2 corresponds to a global spin rotation by 180 
degrees and thus flips back M s , such that, in fact, at the end the order parameter is 
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not affected by O'. As opposed to the honeycomb lattice case, M s changes sign under 
the shift symmetry Di on a bipartite square lattice 18]. In this combined 
shift symmetry D\ leaves the ground state invariant. Since on the square lattice a 
90 degrees rotation O maps sublattices A — )■ A and B — > B, in that case the ground 
state is not affected by a rotation by an angle of 90 degrees. 

Finally, the Hubbard Hamiltonian is invariant under the reflection R at the Xi-axis 
shown in Figure 1. Under this transformation, the fermion operators transform as 



C, 



R^c x R = c Rx . (2.18) 



Since R maps the two sublattices onto themselves, M s remains invariant. 



In [2J, |25j , Yang and Zhang proved the existence of a non-Abelian extension of the 
U(1)q fermion number symmetry in the half-filled Hubbard model. This pseudospin 
symmetry contains U{1)q as a subgroup. The SU{2)q symmetry is realized on the 
square as well as on the honeycomb lattice and is generated by the three operators 



X X 

Q 3 = J2 ^(ir c *t + ^04 - 1) = ^ (2-19) 

X 

The factor (— l) x again distinguishes between the two sublattices A and B of the 
honeycomb lattice. Defining Q 1 and Q 2 through Q ± = Q 1 ± iQ 2 , one readily shows 
that the SU(2)q Lie-algebra [Q a ,Q b ] = i£ a bcQ c , with a,b,c G {1,2,3}, indeed is 
satisfied and that [H,Q] = with Q = {Q l ,Q 2 ,Q 3 ) for the Hubbard Hamiltonian 
with jj, = 0. 

In order to write the Hubbard Hamiltonian Eq. fl2.5p or Eq.f l2.7p in a manifestly 
invariant form under SU(2) S x 577(2) q, we arrange the fermion operators in a 2 x 2 
matrix-valued operator, arriving at the fermion representation 



Cx= ^ V • (2.20) 




The SU (2)q transformation behavior of Eq. ( 12.20]) can now be worked out by applying 
the unitary operator W = exp(id; • Q), 

$C X = W ] C X W = C X Q T , (2.21) 

with 

tt = exp (iu ■ | J G SU{2) Q . (2.22) 

Under an SU(2) S spin rotation, C x transforms exactly like c x , i.e. 

C' x = gC Xl g G SU(2) S . (2.23) 
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Applying an SU(2) S x SU(2)q transformation to Eq. (12.201) then leads to 



®C X = gC x tt T . (2.24) 

Since the SU(2) S spin symmetry acts from the left and the SU(2)q pseudospin sym- 
metry acts from the right onto the fermion operator, it is now obvious that these two 
non-Abelian symmetries commute with each other. Under the discrete symmetries of 
the Hubbard model, C x has the following transformation properties 



A: 


^x 




: 


^x 


= Cox&S, 


O': 


o'rt 

^x 


= {iOl) C xV3 


R : 


^x 


= Crx- 



In terms of Eq. (12.20]) . we are now able to write down the Hubbard Hamiltonian in 
the manifestly SU(2) S , U(1)q, Di, O, O' and R invariant form 

H = -\Y< ^\- c l c v + cjoj + nc\CxCic x \ Trr ;,/>,;. (2.26) 

(xy) x x 

The CT3 Pauli matrix in the chemical potential term prevents the Hubbard Hamiltonian 
from being invariant under SU(2)q away from half-filling. For /j / 0, SU(2)q is 
explicitly broken to its subgroup U(1)q. In addition, the pseudospin symmetry is 
realized in Eq. (|2.26j) only for nearest-neighbor hopping. As soon as next-to-nearest- 
neighbor hopping is included, the SU(2)q invariance gets lost even for \i = 0. The 
continuous SU(2)q symmetry contains a discrete particle-hole symmetry. Although 
this pseudospin symmetry is not present in real materials, it will play an important role 
in the construction of the effective field theory. The identification of the final effective 
fields for holes will lead us to explicitly break the SU(2)q symmetry in Section HI 



2.3 Symmetries of the t-J Model 

Away from half-filling and for U 3> t, the Hubbard model reduces to the t-J model, 
which is defined by the Hamiltonian 

H = vl-t c y + c\c x ) + jJ2S x -S y ~fi 5>+ " 1) \v. (2.27) 

(xy) (xy) x 

Using second order perturbation theory, the antiferromagnetic exchange coupling J 
is related to the parameters of the Hubbard model by J = > 0. Again, t is the 
hopping amplitude, S x is the SU(2) S spin operator on a site x, and fi controls the 
doping with respect to a half-filled system. The projection operator V removes all 
doubly occupied sites from the Hilbert space and hence the t-J model can only be 
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doped with holes. In 26|, the single-hole sector of the t-J model was simulated on the 
honeycomb lattice by using an efficient loop-cluster algorithm. For the construction 
of the effective theory for a hole doped antiferro magnet, the t-J model will serve as 
the microscopic starting point. Except for the SU(2)q symmetry, this model shares 
all symmetries with the more general Hubbard model. 



3 Effective Theory for Magnons 



In this section we investigate the low-energy physics of an undoped quantum antifer- 
romagnet. We will first argue that quantum antiferromagnets are systems featuring 
a spontaneous SU(2) S — » U(l) s symmetry breakdown, which induces two massless 
Goldstone bosons — the magnons. We present the leading-order effective action for 
the pure magnon sector of an antiferromagnet on the honeycomb lattice. In addi- 
tion, a non-linear realization of the spontaneously broken SU(2) S spin symmetry is 
constructed, which will enable us to couple magnons and doped holes in section 4. 



3.1 Low-Energy Effective Action for Magnons 

In quantum antiferromagnets the symmetry group G = SU (2) s of global spin rotations 
is spontaneously broken by the formation of a staggered magnetization. The ground 
state of these systems is invariant only under spin rotations in the subgroup H = 
U(l) s . As a consequence of the spontaneous global symmetry breaking, there are two 
magnons which are described by a unit-vector field 

e(x) = (e 1 (x),e 2 (x),e 3 (x)) E S 2 , e{xf = 1, (3.1) 

in the coset space G/H = SU(2) S /U(1) S = S 2 . Here x = (xi,x 2 ,t) denotes a point in 
Euclidean space-time. The low-energy physics of an undoped antiferromagnet can be 
completely described in terms of the field e(x) which represents the direction of the 
local staggered magnetization. 

Later, we will couple magnons to holes. Since holes have spin 1/2 and are thus 
described by two-component fields, it is convenient to work with a CP(1) representa- 
tion instead of the 0(3) vector representation for the magnon field. We introduce the 
2x2 Hermitean projection matrices P(x) defined by 

P(x) = - [1 + e(x) ■a] = 1 -( } + 6 ^ X \ , 6 ^ ~ ) , (3.2) 

obeying 

P( x )t = p( x ) ) TrP(x) = 1, P{x) 2 = P{x). (3.3) 
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In terms of P(x), to lowest-order in a systematic derivative expansion, the effective 
action for magnons is given by 



S[P] = J d 2 x dt psTrldiPdiP + ^d t Pd t P]. (3.4) 

Here we have introduced two low-energy constants, the spin stiffness p s and the spin- 
wave velocity c. The values of these low-energy constants have been determined very 
precisely using Monte Carlo simulations 27H29|. It should be pointed out that this 
leading-order contribution to the effective action is exactly the same as for an antifer- 
romagnet on a square lattice. Deviations will only show up when higher order terms 
with more derivatives are considered. 

We now discuss how the magnon field P(x) transforms under the various symme- 
tries of the underlying microscopic models. Under global SU(2) S spin transformations 
the staggered magnetization field transforms as 

P(x)' = gP(x)g*. (3.5) 

Note that it is invariant under the Abelian and the non-Abelian fermion number 
symmetries U(1)q and SU(2)q, i.e. 

$P(x) = P(x). (3.6) 

Under the displacement Di and the reflection symmetry R, the sublattices are not 
interchanged such that 

A P(x) = P(x), 

R P(x) = P(Rx). (3.7) 

Under a rotation by 60 degrees, the staggered magnetization vector changes sign, i.e. 
°e(x) = —e(Ox), and therefore 

°P(x) = ^ [1 - e(Ox) • a] = 1 - P(Ox). (3.8) 

Note that in an antiferromagnet on the honeycomb lattice the 60 degrees rotation 
symmetry is spontaneously broken, whereas in an antiferromagnet on the square lat- 
tice, it is the displacement symmetry by one lattice spacing which is spontaneously 
broken. The above transformation property simplifies under the composed symmetry 
O', 

°'P{x) = (ia 2 ) °P{x){ta 2 y = P(Ox)*. (3.9) 

Under time-reversal T, which turns a space-time point x = (xi,X2,t) into Tx = 
(xi,x 2 , —t), the staggered magnetization changes sign and, as a consequence, 

T P(x) = l-P(Tx). (3.10) 
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Since also T is a spontaneously broken symmetry, again it is useful to consider the 
composed transformation T' consisting of a regular time-reversal T and the specific 
spin rotation g = ia 2 - Under the unbroken symmetry T" the magnon field P(x) 
transforms as 

T 'P(x) = (ia 2 ) T P(x)(ia 2 y = P(Tx)*. (3.11) 
The effective action in Eq. fl3.4p is invariant under all these symmetries. 



3.2 Non-linear Realization of the SU(2) S Symmetry 

In order to couple the fermions to the magnons, i.e. to the antiferromagnetic order 
parameter, a non-linear realization of the SU(2) S symmetry has been constructed 
and discussed in detail in 18[. The spin symmetry is implemented on the fermion 
fields by a non-linear local transformation h(x) G U(l) s . This local transformation is 
constructed from the global transformation g G SU(2) S and the magnon field P(x) 
as follows. One first defines a local, unitary transformation u(x) G SU(2) S which 
diagonalizes the staggered magnetization field, i.e. 

u(x)P(x)u(xy = ^(1 + 03) = (J g j , Unix) > 0. (3.12) 

In order to make u(x) uniquely defined, we demand that the element Uu(x) is real 
and non-negative. Using Eq. fl3.2j) and spherical coordinates for e(x), i.e. 

e(x) = (sin^(x) cos<p(x), sin9(x) sin<^(x), cos6*(x)), (3.13) 

one obtains 18l 



u(x) 



1 + ^{x) eAx) - ie 2 (x) 



y/2(l + e 3 (x)) V ~ e i( x ) - ie 2(%) 1 + e 3 (x) 

cos (^r^j sin O^f-J exp(—iip(x)) 

— sin ex P(^V ; '( :E )) cos (^2) 



(3.14) 



Note that the local transformation u(x) rotates an arbitrary staggered magnetization 
field configuration P(x) into the specific constant diagonal field configuration with 
P(x) = |(1 + as). Under a global SU(2) S transformation g the diagonalizing field 
u(x) transforms as 

u(x)' = h{x)u{x)g\ Unix)' > 0, (3.15) 
which implicitly defines the non-linear symmetry transformation 

h( X ) = exp M*) ff3 ) = ( eXP<! ; W) eM _° la{x)) ) e £/(!).. (3.16) 
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The transformation h(x) is uniquely defined since we demand that u n (x)' is again 
real and non-negative. 

The transformation behavior of the field u(x) can be easily worked out from the 
known transformation behavior of P(x). Since u(x) contains only magnon degrees of 
freedom, it transforms trivially under both the Abelian and the non-Abelian fermion 
number symmetries U(1)q and SU(2)q, i.e. 

$u(x)=u(x). (3.17) 
Under the displacement and the reflection symmetry R one finds 

Dt u(x) = u{x), 

R u(x) = u(Rx). (3.18) 

The spontaneous breaking of the 60 degrees rotation symmetry O which takes e(x) 
to —e{Ox) leads to 

°u{x) = t(Ox)u(Ox), (3.19) 



with 

t(x) 



—eAx) + ie 2 {x) 



y/e^x) 2 + e 2 (x) 2 V e i( x ) + ie z( x ) 

— exp(— i(p(x)) 
exp(i(p(x)) 



(3.20) 



Under the combined symmetry O' one finds 

°'u{x) =u{Ox)*. (3.21) 

Since time-reversal T is a spontaneously broken discrete symmetry in an antiferro- 
magnet, it acts on u(x) as 

T u{x) = t(Tx)u(Tx). (3.22) 

On the other hand, the combined time-reversal T' is unbroken and therefore realized 
in a linear manner, i.e. 

T 'u(x) = u{Tx)*. (3.23) 



Finally, we introduce the composite magnon fields v^{x) whose components will 
be used to couple the magnons to the fermions. Using the diagonalizing field u(x), 
we define the composite magnon field 

v^x) = u{x)d^u{x)\ (3.24) 

which under SU{2) S transforms as 

v^x)' = h{x)u{x)g^d ll [gu{x)^h{x)^ = h(x)[v^x) + d^]h(x)l (3.25) 
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Since the field v^(x) is traceless and anti-Hermitean, it can be written as a linear 
combination of the Pauli matrices cr a , 

v^x) = ivfa)<r a , a e {1,2,3}, vfa)EXL (3.26) 

Introducing 



we arrive at 



rf(x)=vl(x)Tivl(x), (3.27) 



Under global SU(2) S transformations the components of transform as 

v l( x Y = V U X ) ~ d M x )i 



v±(x)' = exp(±2ia{x))tf(x), (3.29) 



which indicates that t> 3 behaves like an Abelian U(l) s gauge field, while vt{x) exhibit 
the behavior of vector fields "charged" under U(l) s . The transformation properties 
of the components v^(x) and v^(x) under the discrete symmetries can be worked out 
from the definition of t> M (x) in Eq. (13.241) as well, and are summarized as follows 



>x) + V6d 2 (p{Ux 
)x) + d 2 (f(Ox 



A: 




:*) 


= "fa), 


: 


°vfi 




= \[- v\{Ox) + d lV {Ox) - 




°vli 


» 


= \ [y/3vf(Ox) - V3d 1 <p(Ox 




°vh 




= -vf(Ox) + dMOx), 


0' : 


°'v\i 




= -\[vl{Ox) + V3v 3 2 (Ox)}, 




°'vli 




= \[y/lv\{Ox) -v%(Ox)], 








= -vKOx), 


R : 






= vl(Rx), 




% 3 , 




= -vl(Rx), 








= v*{Rx), 


T : 


T vfi 




= -^ 3 (Tx) + c^(Tx), 




T vh 




= v*(Tx)-d t <p(Tx), 


r ■. 


T 'vh 




= -vKTx), 








= vUTx), 



(3.30) 
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and 



A 
o 



O' 



R : 



T : 



T' : 



Di vt 
o„,± 



o,, ± 



o,, ± 



O' ,± 



o',± 



o'.± 



T ± 
T„.± 



T" ± 
T',± 



-exp(^2icp(Ox))^(vf{Ox) + V3v^(Ox)) : 
exp(T2MOx))l(V3vf(Ox) - v^(Ox)), 

— exp(^2iip(Ox))v f(x), 
-±(v?(Ox) + V3v*(Ox)), 
±(V3vf{Ox)-vT(Ox)), 
-vf{Ox), 

vt(Rx), 

-vf(Rx), 

vt(Rx), 

— Qx^{^2iip{Tx))vf{Tx) 1 
exp(^2iip(Tx))vf (Tx) , 
-vf{Tx), 

vf(Tx). 



(3.31) 



The magnon action of Eq. fl3.4p can now be reformulated in terms of the composite 
magnon field v^(x), 



S\vt 



d 2 x dt 2p s 



v i v i + -^ v t v t 



(3.32) 



At a first glance, the expression v£v~ looks like a mass term of a charged vector field. 
However, since it contains derivatives acting on -u(x), it is just the kinetic term of a 
massless Goldstone boson. 



4 Effective Theory for Magnons and Holes 



In this section we construct a systematic low-energy effective theory for holes coupled 
to magnons. As a first step toward building the effective theory, we identify the cor- 
rect low-energy degrees of freedom that describe the holes. Then the transformation 
behavior of these fermionic fields is investigated in great detail. Finally, the most 
general effective Lagrangian for magnons and holes is constructed. 



4.1 Fermion Fields and their Transformation Properties 

In order to construct the effective theory for hole-doped antiferromagnets, it is es- 
sential to know where the hole pockets are located in momentum space. The disper- 
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Figure 2: The dispersion relation E(k)/t for a single hole in an antiferromagnet on 
the honeycomb lattice simulated in the t-J model for J/t = 2 f^ l. 

sion relation E(k) for a single hole in the t-J model on the honeycomb lattice was 
simulated using an efficient loop-cluster algorithm [26[. The result is shown in Fig- 
ure [21 This simulation clearly shows spherically shaped hole pockets centered around 
(±§^, i^T^) an d (0, ij^) hi the first Brillouin zone. Therefore, doped holes occupy 
the two pockets a and (3 with lattice momenta 

4-7T 

fc« = _fc/9 = ( ,__). (4.1) 
3V3a 

Together with the origin, these two points form a minimal set of three points in 
momentum space. The three points in coordinate space that are related to 0, k a } k? 
by a discrete Fourier transform, define three triangular sublattices A ly A 2 , and A 3 , 
as well as Bi, B 2 , and B 3 on the A- and S-sublattices of the honeycomb lattice. 
The geometry of these six triangular sublattices is illustrated in Figure [3j We now 
introduce fermionic lattice operators with a sublattice index X as an intermediate 
step between the microscopic and the effective fermion fields, 

= u(x)C x , (4.2) 

with x G X, X G {A 1 , A 2 , A 3 , B\, B 2 , B 3 }. The above definition of contains the 
diagonalizing matrix u(x) of Eq. fl3.14p and hence accounts for the non-linearly realized 
SU(2) S symmetry on the effective fermion fields. On even and odd sublattices the 
fermion operator has the following components 




xEX,X e{A 1 ,A 2 ,A 3 }, (4.3) 
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Figure 3: {y4i,y4 2 ,A 3 } and {Bi, B 2 , B 3 } sublattice structure and the corresponding 
primitive lattice vectors. 



and 



x 



U[X) 



-v x 

^1 



xeX,Xe {B 1 ,B 2 ,B 3 }. (4.4) 



Note that with the spontaneously broken spin symmetry only the spin direction rel- 
ative to the local staggered magnetization is still a good quantum number. The 
subscript +(— ) then indicates anti-parallel (parallel) spin alignment with respect to 
the direction of e(x). According to Eqs. (12. 23[) and (I3.15p . under the SU(2) S symmetry 
one obtains 

= u{x)'C' x = h(x)u(x)g^gC x = h(x)*%. (4.5) 



Similarly, under the SU(2)q symmetry one finds 

= $ u (xfc x = u( x )c x n T 



* x n T . 



(4.6) 



The discrete symmetries are implemented on the above fermionic lattice operators 
#5 as 



1 X +Cli 1 



A = 


D t yX 

X 


: 


OyX 

X 


0' : 


O'^X 

X 


R : 


RyX 

X 



t{Ox)9%<t s , 
(ia 2 )^a 3 , 



RX 
Rx ■ 



(4.7) 
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In the effective theory doped holes are described by anticommuting matrix-valued 
Grassmann fields 

consisting of Grassmann field components -?/>± ( x ) instead of lattice operators V^±- We 
also introduce 

= ( 5w ) • <*■*•*>■ <«> 

consisting of the same Grassmann fields as ^ x (:r). Therefore, ^^(x) is not indepen- 
dent of ^ x (x). By postulating that the matrix- valued fields ty x (x) transform exactly 
as the lattice operator one obtains 



SU(2) a 




= h(x)^ x (x), 


SU(2) Q 


^ x (x) 


= ^ x (x)n T , 


A 


D ^ x {x) 


= ^ x (x), 





°m x {x) 


= T(Ox)^ ox {Ox)a 3 , 


0' 


°'* x (x) 


= (ia 2 )V ox (Ox)a 3 , 


R 


R m x {x) 


= V RX (Rx), 


T 


V(x) 


= T(Tx)(ia 2 ) [V x \Tx) T ] a 3 , 




T m x \x) 


= -a 3 [m x {Tx) T ] (<(7 2 )tr(Tx)t 


r 


T 'V x (x) 


= - [^{Txf] a 3 , 




T 'm x \x) 


= a 3 [* x (Txf] . 



Here the transformation behavior under time-reversal T and V is also listed. The 
form of the time-reversal symmetry T for an effective field theory with a non-linearly 
realized SU(2) S symmetry can be deduced from the canonical form of time-reversal in 
the path integral of a non-relativistic theory with a linearly realized spin symmetry. 
The fermion fields in the two formulations just differ by a factor u(x). Note, that 
an upper index T on the left denotes time-reversal, while on the right it denotes 
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transpose. In components the transformation rules take the form 



SU(2) S 




= exp(±ia(x))^p x (x), 


U(1)q 




= exp(iu)^ x (x) , 


A 




= i>± iX (x), 





°^ x (x) 


= ^exp(Tip(Ox))ij° x (Ox 


0' 




= ±4>° x (Ox), 


R 


R ^± (x) 


= ^ x {Rx), 


T 


T ^ x (x) 


= exp(^i(p(Tx))ip±\Tx), 




T ^l\x) 


= — exp(±i(p(Tx))^jj x (Tx), 


r 


T '^(x) 


= -^f(Tx), 




T >f(x) 


= il>£(Tx). 



(4.11) 

Since the spin as well as the staggered magnetization get flipped under time-reversal, 
the projection of one onto the other remains invariant. 

We now want to directly relate the fermion fields to the lattice momenta k a and 
k" , i.e. to the hole pockets a and /3. The new degrees of freedom are thus labeled with 
an additional "flavor" index / G {a,f3}. These fields are defined using the following 
discrete Fourier transformations 



A,/, 



B,f, 



X 



X 



[X 



1 3 

— = ^exp(-ik f v n )^ J 

* n=l 
1 3 

— j= exp(— ik^w n )ip Bn (x), 



(4.12) 



n=l 



where 



Wi 



' 2°' 



V3, 



2 a, 



V5, 



v 2 
w 2 



(a,0), 
(-o,0) 



*'3 
W 3 



2 GE< ) 2 ^ 



(4.13) 



The above vectors connect the discrete three-sublattice structure of A and B in po- 
sition space with lattice momenta k* in momentum space (Figure H]). The fields with 
the pocket (or momentum) index then read 





(x) 


1 

" 71 


exp 


(if) # Al (z) + ^ M {x) + exp (-i^) ^ As (x) 


? 




(x) 


1 

" 71 


exp 


(-if) V Al (x) + ^ Aa (x) + exp (if) ^ M (x) 


•) 




(x) 


i 

" 7! 


exp 


(if) V Bl (x) + ^ 2 (x) + exp (-if) V B3 (x) 


5 




(x) 


i 

" 7! 


exp 


(-if) ty Bl (x) + ty B2 (x) + exp (if) ty B3 (x) 


• (4-14) 
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H> 2 




Vl 



V2 



Figure 4: Sublattice vectors from Eq. \J^.13\j . 



The Fourier transformed matrix- valued fields of Eq. f)4.14p can be written as 



ty AJ (x) = — y^exp(-iA; / t; TO )^ An i 
V 3 n=1 



ty BJ (x) = — ^ exp(-ifc / w n )^ 
a/3 ^"f 

v n=l 

with their conjugated counterparts 

yAf Ux) = ( 

■0_ (x) (x) 



X 



x 



ip+'f(x) — ip B '^ \x) 



ip_{x) tp+ J T (x) 



BJ'h 



(4.15) 



ip B (x) ip_ (x) 
—ip B 'f (x) i/j b '^ (x) 



(4.16) 
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The transformation properties of the fields in Eg. (14.141) are 

h(x)* x ' f (x), 

* x >f(x)n T , 

exp(ik f ai)^ XJ (x), 
exp(-i^)T(Ox)y B 'P(Ox)a 3 , 
exp(i^)r(Ox)y B > a (Ox)a 3 , 
exp{if)T(Ox)^(Ox)a 3 , 
expHf )r(0x)^ A ' a {0x)(T 3 , 
exp(-if)(ia 2 )^(Ox)a 3 , 
exp(*f )(ia 2 )^ B > a (Ox)a 3l 

:2vr' 



SU(2) ■ 

u w \ )s • 




x)' 


OU{Z)q . 




X) 






1 ' 1 

V / 


o ■ 


O^A,a 








T I 




O^B,a 


'x) 






[x) 


O': 


0'^A,a 


[x) 






[x) 






' x) 






[x) 


R : 




[x) 


T : 




[x) 






[x) 


T : 




[x] 






[x) 



,B,a { 

exp(i^)(ia 2 )^(Ox)a 3 , 
exp(-if)(ia 2 )q A > a (Ox)a 3 , 
%? x ' f '(Rx), 

r{Tx){ia 2 ) U Xf '\Tx) T ] a 3 , 



-a 3 



V x ' f '(Tx) T ] (m 2 ) t r(Tx) t , 



CT3 



~^ x f'\T: 
* x >f'{Ta 



^3, 



(4.17) 
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For the Grassmann-valued components we read off 



SU(2) a : 




= exp(±ia(x))ib + (x), 

■L \ \ / / 1 _1_ \ / ' 


U(1)q : 


Q ^ J (x) 


= exp(iu)ip± :l (x) , 


Di : 


Di i>l J {x) 


= expUk-l ' aAib+'Hx) , 


: 


°iPi' a (x) 


= Texp(~i^)exp(TMOx))^(Ox), 




°^(x) 


= Texp{i^)exp{TMOx))^ a {Ox), 




°^(x) 


= Texp{i^-)exp{TMOx))ip^{Ox), 




VJ{x) 


= ^exp(-i^) exp(^fi{p(Ox))-ifj^' a (Ox), 


o ■ 


°'^ a (x) 


= ±exp(-i^U^(Ox), 




°'^{x) 


= ±exp(i^)^' a (Ox), 






= ±exp(i^)ilj^(Ox), 




°'^{x) 


= ±exp(-if)^' a (Ox), 


R : 


R i>i J {x) 


= ipf' f '(Rx), 


T : 




= exp(TMTx))4>± J '\Tx), 






= - exp(±iip(Tx))ip±' f (Tx), 



At the moment, the matrix- valued fermion fields have a well-defined transformation 
property under SU(2)q. Therefore these fields represent both electrons and holes. 
Since we want to construct an effective theory for the t-J model which contains holes 
only, a crucial step is to identify the degrees of freedom that correspond to the holes. In 
order to remove the electron degrees of freedom one has to explicitly break the particle- 
hole SU(2)q symmetry, leaving the ordinary fermion number symmetry U(1)q intact. 
This task can be achieved by constructing all possible fermion mass terms that are 
invariant under the various symmetries. Picking the eigenvectors which correspond 
to the lowest eigenvalues of the mass matrices then allows one to separate electrons 



T' : 




-tf£'\Tx), 
i)l J '{Tx). 



(4.18) 
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from holes. The most general mass terms read 



J2 - Tr[M(* A ' f *<T S * A ' f - * B ' f *<T 3 * B ' f ) + m(^ AJ ^ A ' f a 3 + V B ' f ^ B ' f a 3 )] 



f=a,P 



f=a,P 

o ^ ( ^ 

—M + m 



/=a,/5 



— Al + m 
M + m 



^ A J 

4> BJ 



(4.19) 



The terms proportional to M. are invariant under SU (2)q while the terms proportional 
to m are invariant only under the U(1)q fermion number symmetry. Since these 
matrices are already diagonal, we can directly read off the eigenvalues which are given 
by ±Ai + m. For m = we have a particle-hole symmetric situation. The eigenvalue 
M. corresponds to the rest mass of the electrons, while the rest mass of the holes is 
given by the eigenvalue — M. The masses are shifted to ±M + m when we allow 
the SU(2)q breaking term (m ^ 0), which implies that the particle-hole symmetry 
is destroyed. Hole fields now correspond to the lower eigenvalue — M. + m and are 
identified by the corresponding eigenvectors ^'"(x), ip+'^( x ), ^P- ,a (x), and i^ A '^{x). 
One can show that these hole fields and their conjugated counterparts form a closed 
set under the various symmetry transformations. We can thus simplify the notation, 
since a hole with spin + (— ) is always located on sublattice B (A). Hence, we drop 
the sublattice index and the full set of independent low-energy degrees of freedom 
describing a doped hole in an antiferromagnet on the honeycomb lattice is then given 
by 



ij,«{x)= *{j A > a (x) , if£(x) = ^ p (x) , 



X) 



(4.20) 

Even though SU(2)q will now no longer be considered as a symmetry of the effective 
theory, it was of central importance for the correct identification of the fields for doped 
holes. 
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Under the symmetries of the t-J model the hole fields transform as 



SU{2) S : 


ipl(x)' 


= exp(±ia(x))V'£(a;), 


U(1)q ■■ 


Q ip f ± (x) 


= exp(ia;)'0£(a;), 


Di- 


D ^i(x) 


= exp(i/c^ai)^£(a;), 


Cl : 




= Texp(±i^-TMOx))ij^(Ox), 




V±(x) 


= =Fexp(^i^ ^i(p(Ox))il)°(Ox), 


0' : 




= ±exp(±if )^(Ox), 






= ±exp(^f)^(Ox), 


R : 




= il>i{Rx), 


T : 




= exp(^(Tx))^£ t (Tx), 




T ^± (^) 


= — exp(±i<^(Ta;))^£ (Tx), 


T' : 




= -^(Tx), 




T V£ f (x) 


= ip± (Tx). 



The action to be constructed below must be invariant under all these symmetries. 



4.2 Low-Energy Effective Lagrangian for Magnons and Holes 



The terms in the action can be characterized by the number of fermion fields they 
contain, i.e. 

(4.22) 



^^Lv%V^\ = J d 2 X dt J2 C n,- 



The leading terms in the effective Lagrangian without fermion fields describe the pure 
magnon sector and take the form 



C = 2p s ( vfvt + -^vfv t 



(4.23) 



The leading terms with two fermion fields (containing at most one temporal or two 
spatial derivatives), describing the propagation of holes as well as their couplings to 
magnons, are given by 

f=<x,P 
»=+,- 

+ iX[(A + isa f D 2 )i) f s \v{ + iscr f v s 2 )ip f _ s 
- (vl + isa f v^(D 1 + isa f D 2 )4> f _ s ] 

+ isa f N 2 {^v{v^{ - ^vfc'tf)] . (4.24) 
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Here M is the rest mass and M' is the kinetic mass of a hole, A and K are hole-one- 
magnon couplings, while L, Ni, and N 2 are hole- two- magnon couplings. Note that all 
low-energy constants are real-valued. The sign <jf is + for a and — for /3. We have 
introduced the field strength tensor of the composite Abelian "gauge" field 

ff j {x) = d i v%x)-d j vl{x), (4.25) 

and the covariant derivatives D t and Di acting on i/)±(x) as 

D t i)±{x) = [d t ± ivf(x) - //] ip±(x), 

Diipi(x) = [di ± ivf(x)] ip f ± {x). (4.26) 

The chemical potential \i enters the covariant time- derivative like an imaginary con- 
stant vector potential for the fermion number symmetry U(1)q. It is remarkable that 
the term proportional to A with just a single (uncontracted) spatial derivative sat- 
isfies all symmetries. Due to the small number of derivatives it contains, this term 
dominates the low-energy dynamics of a lightly hole-doped antiferromagnet on the 
honeycomb lattice. Interestingly, for antiferromagnets on the square lattice, a cor- 
responding term, which was first identified by Shraiman and Siggia, is also present 



in the hole-doped case 19]. On the other hand, a similar term is forbidden by sym- 
metry reasons in the electron-doped case [22| . For the honeycomb geometry we even 
identify a second hole-one-magnon coupling, K, whose contribution, however, is sub- 
leading. Interestingly, the field-strength tensor fij appearing in eq. (14.241) and defined 
by eq. (14.251) is not allowed for hole- or electron-doped antiferromagnets on the square 
lattice due to symmetry constraints. 

The dispersion relation for a single free hole of both flavor a and /3 can be derived 
from £2 and is given by 

2M' 

which is just the usual dispersion relation for a free non-relativistic particle. Note that 
V — (piiPz) is defined relative to the center of the hole pockets. Eq. (l4.27p confirms 
that the two pockets a and are of circular shape which is in agreement with the 
result of simulating the one-hole sector of the t-J model on the honeycomb lattice 



E^(p) = M+ 7 ^- + 0(p i ), (4.27) 



26]. 



The leading terms without derivatives and with four fermion fields are given by 

s=+- 

+ GWftftfW + G 3 ^f^W_ s } . (4.28) 

The low-energy four fermion coupling constants G\, G2, and G3 again are real- valued. 
Although potentially invariant under all symmetries, terms with two identical hole 
fields vanish due to the Pauli principle. 
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4.3 Accidental Symmetries 



Interestingly, the leading order terms in the effective Lagrangian for magnons and 
holes constructed above feature two accidental global symmetries. First, we notice 
that for c — > oo and without the term proportional to %K in C 2 , Eq. fl4.23p . Eq. (14.24]) . 
and Eq. fl4.28l) have an accidental Galilean boost symmetry. This symmetry acts on 
the magnon and hole fields as 

G : G P(x) = P(Gx), Gx=(x 1 -v 1 t 1 x 2 -v 2 t,t), 
G ip±(x) = exp(—p{xi + u}ft)if)±(Gx), 
G ip±(x) = exp(p{xi — uH)ip± (Gx), 
G vf(x) = vf(Gx), 
G vl(x) = vl(Gx) - Vivf(Gx), 
G vt(x) = vt(Gx), 

G vf(x) = vf(Gx) - Vivf(Gx), (4.29) 

with 

p{ = M'v l , p( = M'v 2 , ^ = ^. (4.30) 

The Galilean boost velocity v can be derived alternatively by means of the hole dis- 
persion relation in Eq. (14.27]) and is given by Vi = dE^/dp{ for i e {1,2}. Although 
the Galilean boost symmetry is explicitly broken at higher orders of the derivative 
expansion, this symmetry has physical implications, namely the leading one-magnon 
exchange between two holes, to be discussed in the next section, can be investigated 
in their rest frame without loss of generality. 

In addition, we notice an accidental global rotation symmetry 0(7). Except for 
the term proportional to %K , L 2 of Eq. (I4.24p as well as £4 of Eq. (I4.28P are invariant 
under a continuous spatial rotation by an angle 7. The involved fields transform under 
0(7) as 

°^f(x) = exp(isa f ^)i/j f s (0^)x), s = ±, 

°^vi(x) = C0S7 ui (0(7)0:) + SU17 u 2 (0(7)2;), 

°^>v 2 {x) = — sin 7 (0(7)0;) + COS7 v 2 (0(j)x), (4-31) 

with 

0(7)0; = 0(7) (xi, x 2 , t) = (cos7 xi — SU17 x 2 , sin7 x\ + 00S7 x 2 , t). (4.32) 

Here denotes the composite magnon field. This symmetry is not present in the 
A-term of the square lattice. The 0(7) invariance has some interesting implications 
for the spiral phases in a lightly doped antiferromagnet on the honeycomb lattice and 



was investigated in detail in [23 
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Figure 5: Tree-level Feynman diagram for one-magnon exchange between two holes. 

5 One-Magnon Exchange Potentials 

In the effective theory framework, at low energies, holes interact with each other 
via magnon exchange. Since the long-range dynamics is dominated by one-magnon 
exchange, we will calculate the one-magnon exchange potentials between two holes of 
the same flavor a and /3 and of different flavor. 

In order to address the one-magnon physics, we expand in the magnon fluctuations 
nii{x) and 1712(1) around the ordered staggered magnetization 



ex 



f^^lUoW. (5.1) 



For the composite magnon fields this leads to 

[777.2 (x) ± imi(x)\ + O (m 3 ) 



v„ (x) 



2^ 

v l( x ) — 7 — [mi(x)d^m 2 (x) - m 2 (x)d fl m 1 (x)] + O (m 4 ) . (5.2) 

Since vertices with v^(x) involve at least two magnons, one-magnon exchange results 
from vertices with v^(x) only. As a consequence, two holes can exchange a single 
magnon only if they have anti-parallel spins (+ and — ), which are both flipped in 
the magnon-exchange process. We denote the momenta of the incoming and outgoing 
holes by p± and p±, respectively. The momentum carried by the exchanged magnon 
is denoted by q. The incoming and outgoing holes are asymptotic free particles with 
momentum p = (pi,P2) and energy E(p) = M + pf/2M'. One-magnon exchange 
between two holes is associated with the Feynman diagram in Figure El Evaluating 
these Feynman diagrams, in momentum space one arrives at the following potentials 
for various combinations of flavors /, / G {a, (3} and couplings F,Fe {A, K} 

(p{p!]V s F i\p + p_) = V f F l(q) 5{p + +p--pl- PL), F,Fe {A, K} , (5.3) 
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with 



K 2 

v kk(<i) = ~7T- [ 2 (P+i - tofP+z) -Qi + ivfCh] [2(p-i + zoyp-a) + 9i + : 

K 2 2 

^PsQ 



X 



[2(p_i - icrfP-2) + qi- io-fq 2 ] 



V 1k(q) = -9— 2" (9i +^/g2) 2 [2(p_i + za/p_ 2 ) +?i + iff/ft], 

f f Z'Aif r -, 

^AK (?) = l 2 ^ 1 ~ l(T fP-V + ?1 - *°7?2j , 

Vj£a(?) = 7^2~(^ ~ i(T fQ2) 2 [2(p+i - iOfP+2) -qi + icr f q 2 ] , 
iKA 

Vka (q) = -7^- [2(p+i - - gi + *<r/ga] • (5.4) 

We noted earlier that the leading contribution to the low-energy physics comes from 
the A-vertex. From here on, we therefore concentrate on the potential with two A 
vertices only In coordinate space the AA-potentials are given by 

(rlfl\Vll\f + r_) = Vll(f) 6(r + - fl) 5(r_ - t$, (5.5) 

with 

V(L(r) = -§-f ] (r), V[[{f) = exp(2z^). (5.6) 

Here f = f + — r L denotes the distance vector between the two holes and (p is the 
angle between r and the xi-axis. The 5- functions in Eq.f l5.5j) ensure that the holes do 
not change their position during the magnon exchange. It should be noted that the 
one-magnon exchange potentials are instantaneous although magnons travel with the 
finite speed c. Retardation effects occur only at higher orders. 

Interestingly, in the AA channel, one-magnon exchange over long distances between 
two holes can only happen for holes of opposite flavor. For two holes of the same 
flavor, one-magnon exchange acts as a contact interaction. In the next section we will 
concentrate on the long-range physics of weakly bound states of holes and therefore 
we will only consider the binding of holes of different flavor. 



6 Two-Hole Bound States 



We now investigate the Schrodinger equation for the relative motion of two holes with 
flavors a and f3. In the following, we will treat short distance interactions by imposing 
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a hard-core boundary condition on the pair's wave function. Due to the accidental 
Galilean boost invariance, without loss of generality, we can consider the hole pair in 
its rest frame. The total kinetic energy of the two holes is given by 

t= y Tf = y^L = zL. (6.i) 

2M' M' K ' 

f=a,l3 f=a,/3 

We introduce the two probability amplitudes ^ f i(r ) and ^ 2 (f) which represent the two 
flavor-spin combinations and respectively, where we choose the distance 

vector f to point from the /3 to the a hole. Since the holes undergo a spin flip during 
the magnon exchange, the two probability amplitudes are coupled through the magnon 
exchange potentials and the Schrodinger equation describing the relative motion of 
the hole pair is a two-component equation. Using the explicit form of the potentials, 
the relevant Schrodinger equation for two holes reads 

-£A 7 ^exp(-2^) \ / ^(f) \ / Mr) \ f6 2) 

7 ^exp(2^) -^A J ^ * 2 (r) J ~ * \ * 2 (f) ) ' ^ 



with 



A 2 

7 = ■ (6-3) 



2tt P . 

Making the separation ansatz 

*i(r, ip) = Ri(r) exp(imi^), ^ 2 (r, ip) = R 2 (r) exp(im 2 (p), (6.4) 

with r = \f |, and using the Laplace operator in polar coordinates one arrives at the 
coupled equations 

" + \t ~ h ml ) Rl{r) + lM '^ 6XP ( " + mi - m2) ) = M ' ER ^ 

'{^ + \Tr~ ^0 R2{r) + 7M '^ 6XP + mi - m2) ) = M ' ER ^)- 

(6.5) 

The radial and angular part can be separated provided that the condition m 2 -m 1 = 2 
is satisfied. Introducing the parameter m, which is implicitly defined by 

mi = m — 1, m 2 = m + 1, (6.6) 

the radial equations are then given by 

1 ^ D ^ , ~™> R ^ -M'ER^r), 



-(- 


1 d 

r dr 


\dr 2 


} d 2 


1 d 


\ dr 2 


r dr 



- l(m + I) 2 ) R 2 (r) + 7 M'^P = M'ER 2 (r). (6.7) 
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While the other cases would have to be investigated numerically, for m = the two 
radial equations decouple and can be solved analytically. In particular, by taking 
appropriate linear combinations, for m = the two equations can be cast into the 
form 



( ^ + 1 d 1 + (1 + 7M') 1 
\dr 2 r dr J r 2 



( d 2 Id 
\ dr 2 r dr 



+ (1- 7 M')- 



(R^r) + R 2 (r)) = M' E (R^r) + R 2 (r)) , 
(i?i(r) - R 2 (r)) = M'E{R l {r) - R 2 (r)). (6.8) 



Because the two equations are different, but contain the same energy E, one of the 
equations has a vanishing solution. In the first equation the potential always has a 
positive sign and is thus repulsive. In the second equation, on the other hand, the 
potential has a negative sign and is therefore attractive when the low-energy constants 
obey the relation 

M'A 2 

Thus, magnon-mediated forces can lead to bound states only if the low-energy con- 
stant A is larger than the critical value 



A, 



M> 



(6.10) 



Interestingly, the same critical value arises in the investigation of spiral phases in a 



lightly doped antiferromagnet on the honeycomb lattice [23| . There it marks the point 



where spiral phases become energetically favorable compared to the homogeneous 
phase. Here we are interested in the solution of the above system where the first 
equation has a zero solution and the second a non-zero one, i.e. Ri(r) + R 2 {r) = 0. 
Identifying R(r) = Ri(r) — R 2 (r), the second equation takes the form 



d 2 1 d 
dr 2 r dr 



+ (1 - iM') 



1 



R[r) = -M'\E\R(r), 



(6.11) 



where we have set E = —\E\. The same equation occurred in the square lattice case 



19l . |20| and can be solved along the same lines. As it stands, the equation is ill-defined 
because the 1/r 2 potential is too singular at the origin. However, we have not yet 
included the contact interaction proportional to the 4-fermion coupling G3. Here, in 
order to keep the calculation analytically feasible, we model the short-range repulsion 
by a hard core radius r , i.e. we require R(r Q ) = for r < r . Eq.f l6.ll]) is solved by a 
modified Bessel function 



R(r) = AK v (y/M'\E\r). 



v 



(6.12) 



with A being a normalization constant. Demanding that the wave function vanishes 
at the hard core radius gives a quantization condition for the bound state energy. The 
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quantum number n then labels the n-th excited state. For large n, the binding energy 
is given by 

E„~-— ^ V {- m =). (6.13) 



M'r 2 1 V ^HM> - 1 



Like every quantity calculated within the framework of the effective theory, the bind- 
ing energy depends on the values of the low-energy constants. The binding is expo- 
nentially small in n and there are infinitely many bound states. While the highly 
excited states have exponentially small energy, for sufficiently small r the ground 
state could have a small size and be strongly bound. However, as already mentioned, 
for short-distance physics the effective theory should not be trusted quantitatively. 
If the holes were really tightly bound, one could construct an effective theory which 
incorporates them explicitly as relevant low-energy degrees of freedom. As long as 
the binding energy is small compared to the relevant high-energy scales, our result is 
valid and receives only small corrections from higher-order effects such as two-magnon 
exchange. 

Finally, let us discuss the angular part of the wave equation. The ansatz (16.41) 
leads to the following solution for the ground state wave function 

= ( W) ) = BW ( tw ) • «" 4) 

Applying the 60 degrees rotation O and using the transformation rules of Eq. (14.211) 
one obtains 

°y(r,(p) = -V(r,(p). (6.15) 

Interestingly, the wave function for the ground state of two holes of flavors a and j3 
thus exhibits /-wave symmetry^ The corresponding probability distribution depicted 
in Figure El on the other hand, seems to show s-wave symmetry. However, the relevant 
phase information is not visible in this picture, because only the probability density is 
shown. Interestingly, for two-hole bound states on the square lattice, the wave function 
for the the ground state of two holes of flavors a and (3 shows p-wave symmetry, 
while the corresponding probability distribution (whic h ag ain does not contain the 
relevant phase information) resembles d x 2_ y 2 symmetry (l9^ Remarkably, the ground 
state wave function (I6.14p of a bound hole pair on the honeycomb lattice remains 
invariant under the reflection symmetry R, the shift symmetries Di, as well as under 
the accidental continuous rotation symmetry 0(7). 

We would like to emphasize that the /-wave character of the two-hole bound state 
on the honeycomb lattice is an immediate consequence of the systematic effective field 
theory analysis. It seems that the issue of the true symmetry of the pairing state, 



realized in the dehydrated version of Na2CoC>2 x 2/H2O is still controversial [30[. Still, 



it is quite interesting to note that a careful analysis of the available experimental data 



for this compound suggests that the pairing symmetry indeed is /-wave [31 



1 Strictly speaking, the continuum classification scheme of angular momentum eigenstates does 
not apply here, since we are not dealing with a continuous rotation symmetry. 
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Figure 6: Probability distribution for the ground state of two holes of flavors a and (3. 

7 Conclusions 

In complete analogy to our earlier investigations on the square lattice, we have con- 
structed a systematic low-energy effective field theory of magnons and doped holes in 
an antiferromagnet on the honeycomb lattice. Due to the different lattice geometry, 
there are important symmetry differences which have an impact on the allowed terms 
that enter the effective Lagrangian. Interestingly, in contrast to the square lattice 
case, on the honeycomb lattice an accidental continuous spatial rotation invariance 
arises for the leading terms of the low-energy effective Lagrangian. 

As an important result, we have identified the leading magnon-hole vertex which 
results from a term with a single uncontracted spatial derivative. This term, which 
is analogous to the Shraiman-Siggia term on the square lattice, yields a rather strong 
magnon-hole coupling since it appears at a low order in the systematic low-energy 
expansion. As we have investigated earlier, at non-zero hole doping, when A is suf- 
ficiently strong, this term gives rise to spiral phases in the staggered magnetization 
23 . 

In the present work, we have studied the effect of the magnon-hole vertex on two- 
hole bound states. Again in contrast to the square lattice case, it turned out that the 
magnon-hole coupling constant A must exceed a critical value in order to obtain two- 
hole bound states. Our analysis implies that the wave function for the ground state 
of two holes of flavors a and j3 exhibits /-wave symmetry (while the corresponding 
probability distribution seems to suggest s-wave symmetry). This is quite different 
from the square lattice case, where the wave function for the ground state of two holes 
of flavors a and exhibits p-wave symmetry (while the corresponding probability 
distribution resembles d x 2_ y 2 symmetry). 
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We like to stress again that the effective theory provides a theoretical framework 
in which the low-energy dynamics of lightly hole-doped antiferromagnets can be inves- 
tigated in a systematic manner. Once the low-energy parameters have been adjusted 
appropriately by comparison with either experimental data or numerical simulations, 
the resulting physics is completely equivalent to the one of the Hubbard or t-J model. 



Acknowledgments 



F. K. acknowledges that he has made most of his contributions to this paper while 
working in the Condensed Matter Theory Group at the Massachusetts Institute of 
Technology. U.-J. W. likes to thank the members of the Center for Theoretical Physics 
at MIT, where part of this work was done, for their hospitality. C. P. H. thanks the 
Institute for Theoretical Physics at Bern University for their warm hospitality and 
gratefully acknowledges financial support from the Universidad de Colima which made 
a stay at Bern University possible. F.-J. J. is partially supported by NCTS (North) 
and NSC (Grant No. NSC 99-21 12-M003-015-MY3) of R. O. C. This work was also 
supported in part by funds provided by the Schweizerischer Nationalfonds (SNF). In 
particular, F. K. was supported by an SNF young researcher fellowship. The "Albert 
Einstein Center for Fundamental Physics" at Bern University is supported by the 
"Innovations-und Kooperationsprojekt C-13" of the Schweizerische Universitatskon- 
ferenz (SUK/CRUS). 



References 

[1] Z. Y. Meng, T. C. Lang, S. Wessel, F. F. Assaad, and A. Muramatsu, Nature 
464, (2010) 847. 

[2] S. Chakravarty, B. I. Halperin, and D. R. Nelson, Phys. Rev. B39 (1989) 2344. 

[3] H. Neuberger and T. Ziman, Phys. Rev. B39 (1989) 2608. 

[4] D. S. Fisher, Phys. Rev. B39 (1989) 11783. 

[5] P. Hasenfratz and H. Leutwyler, Nucl. Phys. B343 (1990) 241. 

[6] P. Hasenfratz and F. Niedermayer, Phys. Lett. B268 (1991) 231. 

[7] P. Hasenfratz and F. Niedermayer, Z. Phys. B92 (1993) 91. 

[8] A. Chubukov, T. Senthil, and S. Sachdev, Phys. Rev. Lett. 72 (1994) 2089; Nucl. 
Phys. B426 (1994) 601. 

[9] B. I. Shraiman and E. D. Siggia, Phys. Rev. Lett. 60 (1988) 740; Phys. Rev. Lett. 
61 (1988) 467; Phys. Rev. Lett. 62 (1989) 1564; Phys. Rev. B46 (1992) 8305. 



32 



[10] M. Y. Kuchiev and O. P. Sushkov, Physica C218 (1993) 197. 

[11] O. P. Sushkov and V. N. Kotov, Phys. Rev. B70 (2004) 024503. 

[12] V. N. Kotov and O. P. Sushkov, Phys. Rev. B70 (2004) 195105. 

[13] H. Georgi, Weak Interactions and Modern Particle Theory Benjamin-Cummings 
Publishing Company 1984. 

[14] J. Gasser, M. E. Sainio, and A. Svarc, Nucl. Phys. B307 (1988) 779. 

[15] E. Jenkins and A. Manohar, Phys. Lett. B255 (1991) 558. 

[16] V. Bernard, N. Kaiser, J. Kambor, and U.-G. Meissner, Nucl. Phys. B388 (1992) 
315. 

[17] T. Becher and H. Leutwyler, Eur. Phys. J. C9 (1999) 643. 

[18] F. Kampfer, M. Moser, and U.-J. Wiese, Nucl. Phys. B729 (2005) 317. 

[19] C. Briigger, F. Kampfer, M. Moser, M. Pepe, and U.-J. Wiese, Phys. Rev. B74 
(2006) 224432. 

[20] C. Briigger, F. Kampfer, M. Pepe, and U.-J. Wiese, Eur. Phys. J. B53 (2006) 
433. 

[21] C. Briigger, C. P. Hofmann, F. Kampfer, M. Pepe, and U.-J. Wiese, Phys. Rev. 
B75 (2007) 014421. 

[22] C. Briigger, C. P. Hofmann, F. Kampfer, M. Moser, M. Pepe, and U.-J. Wiese, 
Phys. Rev. B75 (2007) 214405. 

[23] F.-J. Jiang, F. Kampfer, C. P. Hofmann, and U.-J. Wiese, Eur. Phys. J. B69 
(2009) 473. 

[24] S. Zhang, Phys. Rev. Lett. 65 (1990) 120. 

[25] C. N. Yang and S. Zhang, Mod. Phys. Lett. B4 (1990) 759. 

[26] F.-J. Jiang, F. Kampfer, M. Nyfeler, and U.-J. Wiese, Phys. Rev. B78 (2008) 
214406. 

[27] U.-J. Wiese and H. P. Ying, Z. Phys. B93 (1994) 147. 

[28] B. B. Beard and U.-J. Wiese, Phys. Rev. Lett 77 (1996) 5130. 

[29] U. Gerber, C. P. Hofmann, F.-J. Jiang, M. Nyfeler, and U.-J. Wiese, J. Stat. 
Mech. (2009) P03021. 

[30] N. B. Ivanova, S. G. Ovchinnikov, M. M. Korshunov, I. M. Eremin, and N. V. 
Kazak, Phys. Usp. 52 (2009) 789. 

[31] I. I. Mazin and M. D. Johannes, Nature Physics 1 (2005) 91. 

33 



